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Abstract 

A novel method of summing asymptotic series is advanced. Such series repeatedly 
arise when employing perturbation theory in powers of a small parameter for compli- 
cated problems of condensed matter physics, statistical physics, and various applied 
problems. The method is based on the self-similar approximation theory involving 
self-similar root approximants. The constructed self-similar continued roots extrap- 
olate asymptotic series to finite values of the expansion parameter. The self-similar 
continued roots contain, as a particular case, continued fractions and Pade approxi- 
mants. A theorem on the convergence of the self-similar continued roots is proved. 
The method is illustrated by several examples from condensed- matter physics. 
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1 Introduction 



The standard difficulty, repeatedly arising in various problems of condensed-matter physics, 
statistical and chemical physics, field theory, and many other problems of theoretical physics 
and applied mathematics, is the necessity of using perturbation theory in powers of a pa- 
rameter assumed to be asymptotically small, while in reality this parameter is either finite 
or even tending to infinity. How then it would be possible to extrapolate the asymptotic 
series of perturbation theory to the finite, or even to infinite values of the parameter? The 
most often used methods of extrapolation are based on Pade approximants [1]. 

But, when extrapolating an asymptotic expansion f\{x) to large values of the parameter 
x, by means of the Pade approximants, one meets the ambiguity, since the straightforward 
use of these approximants yields 

Pm/n(x)~x m - n (x->oo), 

which, depending on the relation between M and N, can tend to either infinity (when 
M > N), to zero (when M < N), or to a constant (if M = N). 

When the character of the large-variable limit is known, one can invoke the two-point 
Pade approximants [1]. However the accuracy of the latter is not high and their definition 
contains several drawbacks, (i) First of all, when constructing these approximants, one often 
obtains spurious poles yielding unphysical singularities, sometimes with a large number 
of poles [1,2]. (ii) Second, there are the cases when Pade approximants are not able to 
sum perturbation series even for small values of an expansion parameter [3]. (iii) Third, 
in the majority of cases, for achieving a reasonable accuracy, one needs to have tens of 
terms in perturbative expansions [1], while interesting problems provide, as a rule, only a 
few terms, (iv) Fourth, defining the two-point Pade approximants, one always confronts 
the ambiguity in distributing the coefficients for deciding which of these must reproduce 
the left-side expansion and which the right-side series. This ambiguity aggravates with 
the increase of the approximants orders, making it difficult to compose two-point Pade 
tables. For the case of a few terms, this ambiguity makes the two-point Pade approximants 
practically unapplicable. For example, it has been shown [4] that, for the same problem, one 
may construct different two-point Pade approximants, all having correct left and right-side 
limits, but differing from each other in the intermediate region by 1000% of uncertainty. 
Hence, in the case of short series the two-point Pade approximants do not allow for getting a 
reliable description, (v) Fifth, the two-point Pade approximants can be used for interpolating 
between two different expansions not always, but only when these two expansions enjoy 
compatible variables [1]. When the expansions have incompatible variables, the two-point 
Pade approximants cannot be defined in principle, (vi) Sixth, interpolating between two 
points, one of which is finite and another is at infinity, one is able to characterize the 
large- variable limit of only rational powers [1]. (vii) Finally, it may happen that in the 
large- variable limit only the power is known, while the amplitude is not. Then the two-point 
Pade approximants cannot be defined. 

There exists a more general approach for extrapolating asymptotic series in powers of 
a small parameter, or a variable, to finite and even infinite values of such variables. This 
approach is based on the self-similar approximation theory [5-17]. In the frame of this 
theory, we have developed the methods of extrapolating asymptotic series by using several 
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types of self-similar approximants, such as optimized approximants, nested exponentials, 
nested roots, iterated roots, and factor approximants [5-17]. 

In the present paper, we advance a novel type of self-similar approximants that may be 
called self-similar continued root approximants, or, for short, self-similar continued roots. In 
a particular case, these continued roots reduce to continued fractions [18] and, respectively, to 
Pade approximants. But, generally, their form is different and not reduceable to continued 
fractions. The self-similar continued roots could be transformed into expressions of the 
type of the numerical nested radicals [19-22], which, however, is not convenient for the 
extrapolation procedure applied to functions. 

In Sec. 2, we explain how the self-similar continued roots arise in the process of the 
self-similar renormalization of asymptotic series and prove the convergence of these root ap- 
proximants. In Sec. 3, we demonstrate, by several examples from condensed-matter physics, 
that the continued roots can be employed as approximants extrapolating asymptotic series 
and providing good accuracy. Possible generalizations for the continued root approximants 
are also mentioned. 



2 Construction and convergence of self-similar 
continued roots 

Assume that we are looking for the solution of a complicated problem characterized by a 
real function f(x) of a real variable x G 1L And let this problem be not solvable explicitly, 
but allowing only for an approximate solution at small values of x, 

f(x) ~ f k (x) (x -+ 0) , (1) 

where it is represented by asymptotic series of orders k = 1, 2, . . ., 

k 

f k (x) = Y,anX n (a = l). (2) 

n=0 

Here, without the loss of generality, we set ao = 1. This is because any function 

g(x) - 9k{x) (x 0) , 
that at small x is given by the series 

k 

9k{x) = g (x)^2a„x n , 

n=0 

can always be written in such a form, where a = 1. Here go(x) is assumed to be known. 
Then the considered function is defined as 

m = 9(x) 



9o(x) 

Respectively, the series at small x take the form 

/ x k 

9k{x_ 
\ 

n=0 



fk{x) = = > a n x , 

9o(x) ^ 
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corresponding to Eq. (2) with do = 1. 

To apply the self-similar renormalization procedure [5-17] to series (2), recall that the 
first-order series, having the linear from f\(x) = l + a±x, transforms into f*(x) = (1 + Aix) s . 
Notice that series (2) can be rewritten in the form 

A(„.(i + ^(i + a,(i + |,(i + ... + ^,(i + ; a.,))))...). 

Then applying the self-similar renormalization sequentially to each of the terms inside the 
brackets, we come to the self-similar continued root 

f* k (x) = (1 + A lX (1 + A 2 x...(l + A k x) s ) s ...) s . (3) 

For instance, in low orders we have 

f*(x) = (1 + A lX ) s , r 2 {x) = (1 + A lX (1 + A 2 x) s ) s , 

fZ(x) = (i + a iX (i + a 2X (i + A 3X yyy . 

Note that, generally speaking, instead of the same power s, we could take different powers 
s k . This, however, would increase the number of parameters that need to be defined, which 
would require to have more information on the sought function. The use of the same powers 
s simplifies the problem. 

One can notice that in the particular case of s = —1, the self-similar continued roots 
reduce to the continued fractions [18] and, respectively, to the Pade approximants. Thus, 

f2k( X ) = Pk/k{x) , 

fik+iix) = Pk/k+i(x) (s = -1) . 

In this way, the self-similar continued roots (3) generalize the notion of continued fractions. 

It is easy to see that if we require the continued roots to be real- valued quantities for all 
x G [0, oo), the parameters A n have to be non-negative, 

A n >0 (n = l,2,...,k). (4) 

The parameters A n are to be defined by the accuracy-through-order procedure, by com- 
paring the like orders in the small-variable expansion for the continued roots (3) with the 
corresponding series (2), 

fi(x)~f k (x) (x^O). (5) 
For instance, in the first order, the expansion of approximant (3) is 

f*(x) ~ 1 + sA lX (x->0), 

which gives 




In the second order, at small x, we have 

sA-t 

f*(x) ~ 1 + sA lX + — (sAx -A ± + 2sA 2 )x 2 , 
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from where we get the same parameter A\ and 

(1 - s)a\ + 2sa 2 



2s 2 ai 



What is left to be defined is the power s. Assume that the large-variable asymptotic 
behaviour of the sought function is known to be 

f(x) ~ Bx p (x ->• oo) . (6) 

In turn, the large-variable behaviour of approximant (3) is 

ftix)*Btffi* (a? -►«>), (7) 

with the amplitude 

B k = f[K n (8) 

n=l 

and the power 

n=l 

Let us assume (which will be proved below) that the approximants converge, under k — > oo, 
so that there exists the limit 

lim (3 k = (3. (10) 

k — >oo 

This requires that the power s, by modulus, be smaller than one, as a result of which, Eqs. 
(9) and (10) yield 

P = ^- s (N<i). (ii) 

Inverting the latter equality gives 

s = T+7 (">-8- 

We have assumed above that the sequence of approximants (3) converges as k increases. 
Now, we provide the proof of this. 

Theorem. Suppose that x lays in a finite interval 

xe[0,L] (L<oo), (13) 
the power s is such that \s\ < 1 and the parameters A n are non- negative and bounded, 

0<A n <M (n = l,2,...). (14) 
Then the sequence of approximants (3) converges, 

im /;(*) = /;,(*). (is) 

fe— too 
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Proof. The self-similar continued root (3) can be reduced to the form of the nested 
radicals 

(1 • (-''2 • (•<•:; • ..-4 SY--Y , 

in which the notation 

n-1 

x n = x^Y[Af n ( Xl = l) 

3=1 

is used, where 

1 - s n_1 

7„ = (n = 2, 3, . . .) . 

1 (1 - s)s n - 1 V ' ' ; 

By the Herschfeld theorem [19], the sequence of the nested radicals, with non-negative x n , 
converges if and only if all x s ™ are bounded. Condition (14) shows that 

n-1 

\\ Af < M 7 " s " , 

3=1 



with 



n 



1 - S 

Taking into account condition (13), we see that the terms x s ™ are bounded for all finite n, 

x s n n < (II) 7 " 8 " . 

When n — >■ oo, under \s\ < 1, then 

<" < (LM) s/(1 - s) (n -)■ oo) . 

Therefore all terms x*™ are bounded for any n, including n — > oo. Hence the sequence of the 
nested radicals converges, as well as the sequence of the self-similar continued roots (3). 

Remark. When considering the large-variable behaviour, it is always possible to treat 
this variable as asymptotically large, but finite, so that condition (13) be valid. 



3 Extrapolation by means of self-similar continued roots 

The self-similar continued roots (3) can be used for approximating functions in the whole 
region of the variable, between zero and infinity. Below, we illustrate that these approximants 
extrapolate, with a good accuracy, the small-variable asymptotic expansion (2) to the large- 
variable region, where the sought function behaves as in Eq. (6). Defining the power s 
according to Eq. (12), we shall calculate the amplitude B k , evaluating the accuracy of the 
found approximations (8) by comparing it with the known value B. 
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3.1 Nonlinear Schrodinger equation 

Many problems in quantum physics, condensed-mater physics, and optics reduce to the 
nonlinear Schrodinger equation of the type 

- \^+(\j + 9M 2 )'l> = Eil>. (16) 

Here we consider the one-dimensional case, where x G (—00, 00). The first term corresponds 
to kinetic-energy operator. The second term represents a harmonic external potential. The 
third term describes the interactions of atoms, typical, e.g., of trapped atomic gases [23-25], 
with a coupling parameter g. 

The ground-state solution represents equilibrium Bose-Einstein condensate. And the 
higher-energy solutions characterize coherent modes [25]. The spectrum of coherent modes 
for Eq. (16) can be written [13] in the form 

K=(n + ^f(a), (17) 

where n — 0, 1, 2, . . . and the function f(a) can be found by means of perturbation theory 
as an expansion 

k 

f k (a) = l + J2a m a m (18) 



m=l 



in powers of the effective coupling parameter 



a _JJr^ (lg) 

1 + 2n y ' K ' 



in which 



1 f°° 

L K(x)e dx ' 



with H n (x) being Hermite polynomial. The first five coefficients in expansion (18) are 

11 13 

ai — 1 , a 2 — - - , a 3 = — , a 4 = - — — , a 5 



32 1 128 1 2048 

At large a, we have 

f(a) ~ ^ a 2/3 (a -)■ 00) . (20) 

Hence, (3 = 2/3, which, in view of Eq. (12), gives s = 2/5. 
The continued roots (3), at large a, lead to the form 

r k {a) ~ B k a^ . (21) 

For the amplitude we find the approximations 



B 2 


= 1.549484 


(3.3%) , 


B 3 


= 1.554034 


(3.6%) , 


B 4 


= 1.539048 


(2.6%) , 


B 5 


= 1.523475 


(1.6%) , 



where in brackets the related percentage errors are shown. 
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3.2 Frohlich optical polaron 

The ground-state energy of the Frohlich optical polaron can be written as a function 

E(a) = -af(a) (22) 

of the effective coupling parameter a. The function f(a) here can be calculated in the 
second-order perturbation theory giving [26] the expression 

f 2 (a) = l + a 1 a + a 2 a 2 , (23) 

with the coefficients 

ai = 1.591962 x 10" 2 , a 2 = 0.806070 x 10~ 3 . 

The strong-coupling asymptotic behavior was found by Miyake [27,28] as 

f{a) ~Ba {a -> oo) , (24) 

with the amplitude B = 0.108513. This gives for the ground-state energy (22) 

E(a) ~ -Be? (a ->■ oo) . (25) 

For function (24), we have /3 — 1, hence s = 0.5. Extrapolating expansion (23) by the 
second-order continued root, we get B 2 = 0.1044, with the error of 3.8%. 

3.3 Fluctuating fluid membrane 

The pressure of a fluctuating fluid membrane, as a function of stiffness g can be represented 
in the form 

Pig) = ^ fig) ■ (26) 

Perturbation theory for the function f(g) has been done [29] up to the six-th order, resulting 
in the expansion 

k 

f k (g) = l + J2 a ng n , (27) 



n=l 



with the coefficients 



a ± = ^ , a 2 = , a 3 = 2.176347 x lO" 3 , 

a 4 = 0.552721 x 10" 4 , a 5 = -0.721482 x 10" 5 , a 6 = -1.777848 x 10~ 6 . 

The pressure of the membrane between hard walls corresponds to the stiffness g — > oo. 
The hard- wall limit was calculated by Monte Carlo techniques [30], yielding 

P(oo) = 0.0798 ± 0.0003 . (28) 

This corresponds to the asymptotic behaviour of the function f(g) as 

f(g) ~ 0.064683# 2 (g ->■ oo) . (29) 
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Therefore here (3 — 2 and s = 2/3. 

Employing the self-similar continued roots, we find the following approximations for the 
hard-wall pressure: 

P 2 *(oo) = 0.047705 (-40%) , 
P 3 *(oo) = 0.061904 (-22%) , 
P 4 *(oo) = 0.072407 (-9.3%) , 
P*(oo) = 0.079569 (-0.29%) , 
P*(oo) = 0.083702 (4.9%) , 
where in brackets the corresponding percentage errors are given. 

3.4 Fluctuating fluid string 

The free energy of a fluctuating fluid string coincides with the ground-state energy of a 
particle in a box [31,32]. The latter, as a function of the wall stiffness g, can be written as 

E(g) = ^ f(g) , (30) 

with the notation 



/G/)=i+£+7\A+£- (3i) 



32 4 V 64 
Perturbative expansion in powers of g yields 



k 



h(g) = l + J2 a n9\ (32) 



n=l 



with the coefficients 

ai = \, a 2 = ±, «3 = ^, a 4 = , 

1 n 1 

Or = , da = , a 7 = , 

5 131072 6 ' 7 16777216 

and so on. 

The case of rigid walls corresponds to g — > oo, which results in 

7T 2 

E(oo) = = 0.077106 . (33) 

128 

This implies the asymptotic behaviour of function (31) as 

f(g) ~ 0.0624998/ (g -> oo) . (34) 
From here, we have (3 = 2 and s = 2/3. 
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We calculate the rigid-wall limit by extrapolating expansions (32) by means of the self- 
similar continued roots. The results for different approximations E£(oo) are presented in the 
Table, together with the related percentage errors 

s S^-EM x 10Q% 



k 


K(oo) 


£k 


2 


0.047705 


-38% 


3 


0.061362 


-20% 


4 


0.070598 


-8.4% 


5 


0.076155 


-1.2% 


6 


0.079097 


2.6% 


7 


0.080336 


4.2% 


8 


0.080533 


4.4% 


9 


0.080141 


3.9% 


10 


0.079540 


3.2% 


11 


0.079199 


2.7% 


12 


0.079123 


2.6% 


13 


0.078363 


1.6% 



Table: Self-similar continued root approximants of different orders k for the hard-wall energy 
El(oo), with the related percentage errors. 

3.5 Discussion and possible generalization 

The problem is considered of extrapolating asymptotic series in powers of a small variable 
to large values of this variable. A new structure is shown to arise as a result of the self- 
similar renormalization, giving the self-similar continued root approximants. The use of 
these approximants for extrapolation procedure requires the knowledge of the power for the 
large- variable behaviour. The convergence of the sequence of these approximants is proved. 
Several examples from quantum condensed matter physics illustrate good accuracy of the 
approximations. These examples show that convergence is not necessarily monotonic, but 
can be oscillating. Nevertheless, the convergence theorem tells us that the sequence of these 
approximants does converge. 

The method is applicable, when all parameters A n in form (3) are non-negative. It 
may happen that for particular cases some of these parameters occur to become negative, 
which would yield to complex expressions for the roots. In that situation, the method 
cannot be used directly, but one can proceed by taking the highest available real-valued root 
approximation and using the following terms for constructing corrected approximants as is 
explained in Ref. [17]. 

We have checked the applicability of the method for a number of other problems. As a 
rule, the low-order approximants are usually real, giving good accuracy, but may become 
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complex for higher-order approximants. For instance, if we try to approximate, by continued 
roots, the statistical sum of the zero-dimensional oscillator as a function of a coupling pa- 
rameter [33], when (5 = —1/4 and s = —1/3, we get in the second order the strong-coupling 
amplitude 0.970, whose error is 5%. The third and fourth approximants are real, though 
slightly worse than the second one. However, the fifth approximant is complex. 

In another example, we calculate the expansion factor of a three-dimensional polymer 
chain as a function of an effective coupling parameter [34,35], when (3 = 0.3544 and s = 
0.261666. The second-order approximant for the amplitude of the strong-coupling behavior 
gives 1.554, which is within an error of 1.5%. The third-order approximant is yet real, but 
the fourth is complex. Thus, in the above two examples, we need to limit calculations by 
the low-order approximants. However, it is worth recalling that the majority of the most 
interesting problems in condensed-matter physics rarely enjoy the luxury of having many 
expansion terms. The standard situation is when one is able to derive just a few lowest 
terms of perturbation theory. 

In the main expression (3) that has been used throughout the paper, we have taken the 
same powers s. Strictly speaking, it has been admissible to write a more general expression 

f* k (x) = (1 + A lX (1 + A 2 x...(l + A k x) Sk ) Sk -' . . .) S1 , 

with different powers s n . But then, it would be necessary to have more information on the 
large-variable behaviour in order to define all these different powers. 

Concluding, we have suggested a novel method for extrapolating the small-variable 
asymptotic series to the large values of the variable. The method produces a convergent 
sequence of approximants and can be used for extending the validity of perturbation theory. 
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